Continuum Description of the Cytoskeleton: Ring Formation in the Cell Cortex 
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Motivated by the formation of ring-like filament structures in the cortex of plant and animal cells, 
we study the dynamics of a two-dimensional layer of cytoskeletal filaments and motor proteins near a 
surface by a general continuum theory. As a result of active processes, dynamic patterns of filament 
orientation and density emerge via instabilities. We show that self-organization phenomena can lead 
to the formation of stationary and oscillating rings. We present state diagrams which reveal a rich 
scenario of asymptotic behaviors and discuss the role of boundary conditions. 
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The cytoskeleton, an organized network of filamentous 
proteins, is an essential component of all eukaryotic cells. 
It plays a major role in morphogenesis, transport, motil- 
ity and cell division The key components of the cy- 
toskeleton are actin filaments and microtubules which are 
long elastic protein filaments that reach lengths of sev- 
eral /^m. They have chemically and dynamically distinct 
ends and hence are structurally polar. Motor proteins 
consume chemical fuel to generate directed motion on 
these filaments and hence can induce stresses in the cy- 
toskeletal network. 

An important cytoskeletal structure in various cell 
types are cortical rings formed by bundles of filaments 
which wrap around the cell. Such rings form within the 
cell cortex, a thin layer of filament network located close 
to the cell membrane. In animal cells a contractile ring 
containing actin filaments, cleaves the dividing cell by 
generating a constriction [lj. In higher plant cells, mi- 
crotubules form cortical arrays (CA) between consecu- 
tive cell divisions [2|, with a preferred orientation in the 
azimuthal direction, see Fig.^a). Shortly before enter- 
ing division, the CA is dynamically reorganized and the 
preprophase band (PPB) is formed 0]. It consists of a 
ring-like bundle of microtubules that determines the lo- 
cation and orientation of the future division plane of the 
cell 0, see Fig. ^b). In addition to single rings, double 
rings have been observed in certain cells p|. 

It has been experimentally observed in vitro that 
mixtures of microtubules and motor proteins can self- 
organize into a number of patterns, such as bundles, 
asters and vortices |E |E] . Similar patterns where found 
in molecular dynamics simulations [E Q. Theoretical 
descriptions of cytoskeletal dynamics are either based 
on microscopic models for the motor-filament interac- 
tions n n m ei or consider the behavior on large 

length and time scales in terms of a continuum the- 
ory [H EE EE EE EE E3 • In contrast to microscopic de- 
scriptions, which rely on specific and potentially incom- 
plete microscopic pictures of the nonequilibrium driving 
forces, a phenomenological description captures generic 
conditions for the formation of patterns through insta- 



FIG. 1: (color online) Density and orientation of cortical fila- 
ments. (a,b) Fluorescence microscopy images of microtubule 
patterns in Tobacco BY2 suspension cells, (a) Between cell di- 
visions microtubules form a cortical array and orient preferen- 
tially perpendicular to the cell's long axis, (b) Shortly before 
cell division microtubules condense into a ring-like structure, 
the Pre-Prophase Band (PPB). (c) Representation of a sta- 
tionary solution to the dynamic equations resembling a 
PPB. The density of filaments is represented by the grey level 
on the cylinder. The blue bars indicate the nematic order of 
filaments by their orientation and length. (a,b) courtesy of 
Jan Vos, Wageningen University. Scale bars are 10 /im. 



bilities. In the present work we develop a general coarse 
grained description of the two-diemensional cell cortex 
based on symmetry arguments and extend previous work 
to include nematic order. 

We focus on the dynamical behaviour of a cortical 
layer of filaments at scales large compared to the fil- 
ament lengths and the layer thickness, where a coarse 
grained, two-dimensional description is suitable. We as- 
sume that the state of the system is essentially specified 
by the scalar mass density c(r, t), the polarization p(r, t), 
a vector along the average direction of the filaments, and 
q(r,t), a second rank traceless symmetric tensor that 
measures orientational order. These fields represent slow 
hydrodynamic modes as they are linked either to a con- 
servation law (the density c) or to a spontaneously broken 
symmetry (the fields p and q). Other degrees of freedom 
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in the system are assumed to relax more rapidly and are 
thus ignored [lif . 

The choice of the appropriate dynamic equations for 
these fields is governed by the following principles (i) all 
dynamics is overdamped as its origin lies in the motion of 
colloidal sized particles in a viscous background fluid, (ii) 
all mass transport is due to stresses generated in the sys- 
tem, and (iii) all driving terms are allowed that respect 
the rotational invariance of space. This is consistent with 
the fact that they are caused by local interactions be- 
tween the microscopic components. For the mass density 
we have the continuity equation 



d t c + V • j = s 



(1) 



where j is the filament current and s represents source 
and sink terms describing the polymerization and depoly- 
merization of filaments. The filament current is a conse- 
quence of inhomogeneities in the mechanical stress in the 
filament system. For a layer of cortical filaments close to 
the cell surface moving against an immobile background 
fluid and assuming local isotropic friction, we thus write 



(2) 



where a is the stress tensor and r\ an effective friction 
coefficient. The dynamic equations for the non-conserved 
fields p and q can be expressed as d t p = t/>, and d t q = 
where ^ is a vector and a; is a symmetric traceless tensor. 

These equations are closed by assuming that the evo- 
lution of the system is completely determined by the 
present state of the system. This implies that the stress 
tensor <x and the "velocities" of the polarization and the 
nematic order parameter xji and uj can be expressed in 
terms of the dynamic fields c, p and q. We then ex- 
pand these driving forces in terms of perturbations of 
the dynamic fields and their spatial derivatives around a 
reference state that we assume to be homogenous. 

We describe the cell cortex as a two-dimensional active 
gel, embedded in three dimensional space along the cell 
surface. Motivated by the morphology of plant cells, we 
wrap the 2d gel on a cylinder. Experimental evidence 
indicates that in the CA filaments orient in a non-polar 
way 0. It has been suggested that this holds for the 
PPB as well Therefore, we neglect the polariza- 

tion p and only keep the nematic order q. In situations 
where the total mass of polymerized microtubules is con- 
served (such as the transition from CA to PPB j!]]), the 
source term s can be neglected. For simplicity, we focus 
on states of rotational symmetry around the cylinder axis 
and we neglect effects of the surface curvature. 

The cylinder axis is the x-axis of our coordinate frame, 
the second dimension described by the y-eods is wrapped 
around the cylinder with radius R, see Fig.^c). The dy- 
namics of configurations which are rotationally invariant, 
obey a set of equations projected on the cylinder axis, 



which can be derived by symmetry arguments \l2\ Il5j : 

7]- l u xx = A lC + A 2 c 2 + A 3 cdlc + A 4 (d x cf + A 5 qJ£) 
U xx = -Biq xx + B 2 d 2 x q xx + B 3 d x q xx + B 4 d 2 x c (4) 

„2 i „2 



+B 5 (q xx + q xy )qxx 

—B\q xy + B 2 d 2 q xy + B 3 d x q xy 
+ B 5(ql x + qL)qx V 



(5) 



where we have neglected higher order terms. It has been 
shown that the terms in Eq. 10-lJHJ) arise from coarse- 
grained microscopic models and can be related to fila- 
ment sliding induced by mot or p roteins in an immobile 
viscous background [jl lrilll5ll27l |. An example is the con- 
tractile tension given by the term A 2 c 2 which has been 
discussed in 

Dimensional analysis reveals that these equations con- 
tain six independent parameters. The values of only three 
of these turn out to influence the dynamics qualitatively; 
they can be identified with A 2 , B 2l and F — —Ac,B±. 
The other parameters describe nonlinear terms which 
are required to stabilize the dynamics. Stability requires 
A 3 > 0, A 4 < 0, B 3 < 0, and B 5 < 0. The specific 
values of these parameters do not qualitatively influence 
the asymptotic dynamics. 

In order to get some general insight into this problem, 
we first consider the somewhat artificial case of periodic 
boundary conditions in x-direction with period L. In 
this case, the homogenous isotropic state, i.e., c(x) = cq, 
qxx(x) = q X y{x) = for all x, is a stationary solution of 
the dynamic equations ©-(JSJ. We first perform a linear 
stability analysis of this state. The linearized equation 
for the perturbation Sq xy decouples from the correspond- 
ing equations for the perturbations Sc and Sq xx . The field 
q xy becomes unstable for B\ < 0. This instability corre- 
sponds to the isotropic-nematic transition. Here, we fo- 
cus on the instabilities of the isotropic homogenous state 
with Bx > 0. 

The linearized dynamics for Sc — ^2 n Cke 27Tmx ^ L and 
hxx = Y, n qxx,ke 2 ™ lx/L is given by 



d_ ( c fc 

dt \ q xx ,k 

-Dk 2 - A 3 c fc 4 
-B 4 k 2 /2 



(6) 



+A 5 k 2 
-Bx - B 2 k 2 4 



B 3 k A 



('k 

qxx,k 



where D = —{A\ + 2cqA 2 ) is an effective diffusion con- 
stant. Figure |21 represents the region of stability of the 
homogenous isotropic state as a dark grey area for suf- 
ficiently large values of the parameters D and B 2 . The 
parameter —D is related to the contractile tension in the 
system. For sufficiently small D, the system tends to un- 
dergo density instabilities which lead to ring solutions. 
The parameter B 2 characterizes a length scale related to 
the orientational order. For sufficiently small B 2 orien- 
tational instabilities occur. The corresponding unstable 
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FIG. 2: Schematic state diagrams of an active gel layer in 
cylindrical geometry with periodic boundary conditions in the 
^-direction, are shown for two values of _F as a function of two 
dimensionless parameters B2 and D. Regions of linear stabil- 
ity of homogeneous filament distributions are shaded in dark 
grey. Outside these regions, the topology of the state dia- 
gram as determined numerically is indicated by dotted lines. 
Asymptotic states include chevron patterns and single or mul- 
tiple rings that are either stationary or do oscillate. Regions 
of oscillating solutions are shaded in light grey. The long 
dashed line indicates a Hopf bifurcation determined by linear 
stability analysis. Examples for stationary chevron and ring 
patterns are displayed in Fig. [3] Parameters are defined as 
D = -(Ai +2coA 2 ), B 2 = B 2 /{B 1 L 2 ), F = F/(B%L 4 ) and 
A'2 — A2/(BiL 4 ). Parameters values are en = coL = 0.5, 
L(-2B 1 /A 1 ) 1 / 2 =5,Ai = A 1 /(B 1 L 2 ) = -2, ^/(BiL 6 ) = 
0.1, Ai/iBiL 6 ) = -1, A 5 /B 1 = 1, S 3 /(SiL 2 ) = -0.05, 
B5/B1 = -10 and F = -1 for (a) and F = I for (b). L 
is the period of the system. 



mode consists of periodically alternating regions with ne- 
matic order at an angle of 7r/4 and — 7r/4 with the re-axis. 
These structures are reminiscent of chevrons in nematic 
liquid crystals, cf. Fig. G^a) For F < 0, the ho- 

mogeneous state becomes unstable by either stationary 
instabilities or by Hopf bifurcations which lead to trav- 
eling wave solutions. Because of symmetry, these waves 
can propagate along the x-axis in both directions. Hopf 
bifurcations occur if (D - B 2 ) 2 + 4F < 0. For F > in- 
stabilities towards inhomogeneous states are always sta- 
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FIG. 3: (color online) Examples of stationary solutions to 
the dynamic equations with periodic boundary conditions, (a) 
Chevron pattern with homogeneous filament density for B 2 = 
-1.25, D = 0.75, and F = -1. (b) Multiple ring pattern for 
B 2 = -1.5, D = 1.5, and F = -1. 



tionary. The instabilities towards stationary states lead 
to periodic patterns of filament accumulation combined 
with nematic order, parallel or perpendicular to the x- 
axis, at the maxima of the filament density, see Fig.^b). 
For decreasing B2 the number of rings increases as the 
corresponding length scale decreases. In order to obtain 
the full phase diagrams displayed in Fig. 2, we solved 
the dynamic equations (3)-(5) numerically with periodic 
boundary conditions 23] . 

Periodic boundary conditions are inappropriate to de- 
scribe cells. An important boundary condition is the zero 
flux condition j'p^l = 0, which implies that the cortical 
material cannot leave the cell. Additional boundary con- 
ditions have to be specified to fully define the asymp- 
totic solutions. It is observed that the end-faces of inter- 
phase plant cells are essentially free of microtubules. This 
is consistent with the absence of imposed order at the 
boundaries and we therefor set qo,L — 0. In the follow- 
ing we use for simplicity the boundary conditions j = 0, 
q = and d x c = d x q xx = d x q xy = at both x = and 
x = L. 

Figure njc) shows a steady state solution of Eqs. (J3J- 
|JSJ with zero flux boundary conditions. This solution 
corresponds to a stationary filament ring with filaments 
that are oriented in y-direction and localized in a ring-like 
pattern wound around the cylinder, strongly reminiscent 
of the PPB. The corresponding stress is anisotropic and 



for Oyy > a xx the ring solution is contractile. Additional 
asymptotic solutions are similar to the solutions shown in 
Fig. |3 with periodic boundary conditions. Furthermore, 
oscillatory solutions appear along a line of Hopf- bifurca- 
tions. An example of such a state with no flux boundary 
conditions is shown in Fig. ^ The oscillating solutions 
for no flux boundary conditions correspond to traveling 
waves in the case of periodic boundary conditions. Ad- 
ditional states can be generated via the boundaries by 
imposing other boundary conditions. 

In conclusion we have presented a generic approach 
to the dynamics of a cortical layer of an active gel of 
cytoskclctal filaments near the cell surface. We have 
focussed on a situation which can describe the forma- 
tion of filament rings such as the PPB in plant cells and 
the contractile ring in animal cells. Ringlike cytoskele- 
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FIG. 4: (color online) Snapshots at different times of an os- 
cillating filament ring for no-flux boundary conditions and 
Bi — 1, D — 0.25 and F = 1. The pattern consists of 
a filament ring which forms near the center of the cylinder 
(a) and moves towards the right pole (b,c), where the ring 
disappears. Simultaneously, a new ring is formed near the 
center (d) which subsequently moves to the opposite pole. 
The whole process is repeated periodically. The times cor- 
respond to phases <j> = cut of the oscillation with </> — (a), 
<f> = (4/36)2tt (b), <f> = (7/36)2tt (c) and <j> = (15/36)27r(d). 



tal patterns occur also in prokaryotes. An example is 
the Z-ring formed by the tubulin analog FtsZ |24|. We 
have analyzed the formation of patterns via dynamic in- 
stabilities resulting from active processes such as the ac- 
tion of motor proteins and the polymerization and de- 
polymerization of filaments. Of particular interest is the 
spontaneous formation of filament rings which provides a 
possible mechanism for the formation of the PPB. Such 
ring structures can be contractile and could induce con- 
striction of a deformable cylinder. Furthermore, we have 
found solutions corresponding to filament rings which 
move periodically towards both ends of the cylinder. Os- 
cillating actin rings have been observed in lymphoblasts 
treated with nocodazole to remove the microtubules |25| 
and more recently in fragments of fibroblasts |2f3 | . Quali- 
tatively, they exhibit the same type of oscillatory pattern 
as described in Fig. 0] 

Our description assumes that filaments are short as 
compared to the cylinder radius R. In the PPB, micro- 
tubules could have a length of several tens of /xm. In 
this case, the bending elasticity can become important 
and could lead to anisotropies in our description since az- 
imuthal orientation would be energetically disfavored. In 
order to connect our macroscopic approach with molecu- 
lar properties of filaments and associated proteins, more 
microscopic descriptions are valuable. For example, neg- 
ative values of the parameter D in the phase diagram of 
Fig. 2 can result from contractile stresses caused by the 
action of motor-aggregates 0,0. Similarly, the param- 
eter B2 characterizes a length scale associated with ori- 
entational order, which depends for example on filament 
lengths and packing density. Such relations between 
microscopic models and phenomenological descriptions 
could help to identify microscopic mechanisms, which un- 
derlie the formation of cortical rings in cells. 
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